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Abstract 



■ Recent advances in the photon tagging facilities together with the novel, high resolution fast 

calorimetry made possible to perform photoproduction cross section measurements of pseudoscalar 
mesons on nuclei with a percent level accuracy. The extraction of the radiative decay widths, needed 



for testing the symmetry breaking effects in QCD, from these measurements at small angles is done 
by the Primakoff method. This method requires theoretical treatment of all processes participating 



> 

| in these reactions at the same percent level. The most updated description of general processes, 



including the nuclear coherent amplitude, is done in our previous paper. In this work, based 
on the framework of Glauber multiple scattering theory, we obtain analytical expressions for the 
incoherent cross section of the photoproduction of pseudoscalar mesons off nuclei accounting for the 
mesons' absorption in nuclei and Pauli suppression at forward production angles. As illustrations 
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of the obtained formulas, we calculate the incoherent cross section for photoproduction from a 
closed shell nucleus, 16 0, and from an unclosed shell nucleus, 12 C. These calculations allow one 
to compare different approaches and estimate their impact on the incoherent cross section of the 
processes under consideration. 
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I. INTRODUCTION 



Recent years have seen a renewed interest in the study of photoproduction of pseudoscalar 
mesons off complex nuclei at few GeV energies due to the advent of high duty-cycle electron 
accelerators together with new photon tagging facilities and novel experimental technolo- 
gies. It has now become possible to perform photoproduction cross section measurements 
of pseudoscalar mesons off nuclei with a percent level accuracy. 

Photoproduction of pseudoscalar mesons off nuclei at small angles, 

1 + A^Ps + A (Ps = ir°,ri,Ti'), (1) 

is provided mainly by the possibility of extraction of their radiative decay widths from 
experimental data. These reactions at small angles proceed via production in the Coulomb 
field of a heavy nucleus, the Primakoff effect [lj. On the other hand, photoproduction in 
the nuclear strong field complicates the extraction of decay widths from experimental data 
and requires detail theoretical examination. 

The precision measurement of 7r° radiative decay width via the Primakoff effect in the 
photoproduction of pseudoscalar mesons off nuclear targets has recently been completed at 
Jefferson Lab The energy upgrade of the CEBAF accelerator to 12 GeV will provide a 
great opportunity to extend the current program to the r\ and rj' mesons. 

In a recent article {3} we discussed that a precise knowledge of the radiative decay widths 
of light pseudoscalar mesons provide an accurate test of chiral anomalies and mixing effects 



that are due to isospin breaking by the difference of the masses of light quarks, 
great importance to understand fundamental properties of QCD at low energies 



lis is of 
(J. The 



precision extraction of the decay widths requires state-of-the-art theoretical descriptions of 
all processes participating in the reaction. In our recent work [3!] we considered explicitly the 
photoproduction of pseudoscalar mesons in the Coulomb field of nuclei and their coherent 
photoproduction off the nuclear strong field. In the present work we will consider the 
incoherent photoproduction of pseudoscalar mesons off nuclei, i.e., 

7 + Ps + A'. (2) 

Here A' is the target final state including all possible nuclear excitations and its breakup. 

fl 

As in our previous work |3j , here we will discuss the incoherent photoproduction of ir° meson 
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bearing in mind that the obtained expressions, after obvious modifications, can be applied 
to photoproduction of the other pseudoscalar mesons. 

The necessity for such investigations is two-fold. As mentioned above, the advent of 
precise experimental data demands the accuracy of theoretical predictions to be at relevant 
levels. Thus, the widely used expressions for the incoherent photoproduction of pseudoscalar 
mesons j3, S] cannot be considered as satisfactory. The reason for this is that the nucleon 
spin-nonflip amplitude 1 is zero at zero production angle due to conservation of angular 
momentum. 



As was first shown by Faldt [U|], such behavior of the nucleon amplitude requires partic- 
ular consideration for the coherent photoproduction process off nuclei. On the other hand, 
in the incoherent production, as we will show later, similar corrections to the cross section 
take place and their neglect leads to meaningless results for the incoherent photoproduction 
at forward angles. 

Despite the fact that the incoherent production off nuclei in non-diffractive processes 



was considered in the works of Faldt 
since the expressions obtained in Refs. 



13|. we are compelled to study this process again 
131 ] are not suitable for the production at small 
angles 2 . Pauli suppression is at work for production at small angles and its correct accounting 
is crucial for the process under consideration. 



II. THE INCOHERENT PHOTOPRODUCTION OFF NUCLEI 

In the Glauber multiple scattering theory [14J] the amplitude of the process as a result 
of which the nuclear wave function in the initial state <&i(ri, ^...ta) transforms to the wave 

1 The photoproduction of pseudoscalar mesons off the nucleon is described by four independent ampli- 
tudes [9(. In the coherent photoproduction off spinless nuclei only the spin-nonflip amplitude survives. As 
for the incoherent photoproduction, this amplitude is the dominant one [lo| , though all four amplitudes 
are involved in the process. 

2 For instance, the incoherent cross section would be strongly suppressed as a result of Pauli exclusion 
principle, which is not the case in Ref. [l3| (see, for instance, expression (33b)). 
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function of the final state $/(ri, r2...r A ) is 

F fl (q,A) = % ^-J d 2 be i ^d 3 r 1 ...d 3 r A ^(f 1 ,f 2 ...f A ) 

A 

y~] Tj (b, ft . . .r A ) ®i(rx,r 2 ...r A ), 



A 

X 

3=1 



r j (b,f 1 ...f A ) = Y p {b-Sj)e l 

A 

x JJ[l-r s (6-5i)^-^)]- (3) 

Here the two-dimensional vectors b and s*j are the impact parameter and the nucleon's 
transverse coordinate, respectively; z\ is the nucleon's longitudinal coordinate in the nucleus; 
q and A are the transverse and longitudinal components of the transferred momenta: q 2 = 
4kpsm 2 (9/2), A = M 2 /(2k). Here k and p are the photon and meson momenta and M is the 
meson mass. The profile functions T PtS (b— s), by definition, are the two-dimensional Fourier 
transformations of the elementary amplitudes for the photoproduction of pseudoscalar meson 
off the nucleon, f p = f(^ + N — ► Ps + N), and elastic meson-nucleon scattering, f s = 
f(Ps + N -> Ps + N): 

r PlS (b- = ^ / fpA*)** ( 4 ) 

The summed cross section, including both the coherent and the incoherent production pro- 
cesses, can be obtained if one uses the closure approximation, i.e., assumes that in the limits 
of the produced particle energies (usually 50-100 MeV) the final states of the nucleus forms 
a complete set: 



= (5) 



Under this assumption the summed cross section depends only on the ground state wave 
function Q: 

— = VlF- f | 2 



dtt 

f 

k 2 



e i ^- i ' ) d 2 bd 2 b'd 3 r 1 ...d 3 r A \^ i (f 1 ...f A Y 2 



(27T) 2 

x Y, T i&n.--7A^)T*,(b',r 1 ...r A ,A). (6) 



3,3' 
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Let us assume that the ground state of the nucleus can be described by means of the 
independent particle model. Introducing the single particle density p(r) we have 

\Mri,...r A )\ 2 =11^; [p(rld*r = A. (7) 



k=l 



Separating the contributions of diagonal and non-diagonal terms in the sum in Eq. [6l one 



gets 
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15|: 



do a 



e i *$^d 2 bd 2 b'd 2 sdz 



dn (2tt) 2 

x r p (b - s)r;(b' - s) P (s, z)e~ E ^\ 

POO 

E(b,b',z) = / [T s (b-s) + T:(b'-s) 

J z 

- T 8 (b- s)T* s (b' - s)]p(s,z')d 2 sdz', 
dond (A -I) k 2 

dQ A (2tt) 2 

X 23? J e i ?-(P-b')+iA(z 2 -z 1 ) 

x 9(z2 — zi)d 2 bd 2 b' d 2 s\d 2 $2dz\dz2 

x T p {b-s x )Y* p {b' -s 2 )[l-Y s {b' -s 2 )] 

x p(s 1 ,z 1 )p(s 2 ,z 2 )e- E ^'' z ^\ 

E(b,b', zi, z 2 ) = / T s (b — s)p(s, z)d 2 sdz 

J Zl 

POO 

+ / [r s (b-s) + r:(b'-s) 

J Z2 

- T s {b - s)F* s (b' - s)]p{s,z)d 2 sdz. (8) 

These expressions are general and can be applied to diffractive production off nuclei (for 
instance, in photoproduction of vector mesons V(p, u, <p, ip) ) and to non-diffractive processes 
like production of pseudoscalar mesons. 

The incoherent processes at large momentum transfer were considered in detail in 



Ref. 



1 2t] 3 . Our main goal here is to investigate the incoherent photoproduction of pseu- 



doscalar mesons at small momentum transfers. In this kinematic region one can neglect the 



In Rcfs. the final formulas have been cited for the production processes with the same absorption 

in nuclei for the initial and final hadrons, which significantly facilitates the calculations, but it is not 
appropriate for photoproduction. 
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mesons' multiple scattering to non-forward angles, the effect important at high momentum 
transfers [16]. In this approximation the differential cross section Eq. [8] takes the form: 

da da d da nd 
dVt ~ dtt dtt 



+ 



p(s, z)\h(s, z)\ 2 d 2 sdz 
(A- I 



A 



h(s, z) 



ik 



2tt 
x exp 



p(s, z)h(s, z)d 2 sdz 
e l ^ b+iAz d 2 bT p (b - s) 

T s (b — s')p(s , z')d 2 s'dz' 



(9) 



The second term in this expression includes the coherent cross section. To obtain the 
incoherent cross section one has to subtract it from the summed cross section (Eq. [9]): 

define 



dQ 



1 

A 



p(s, z)\h(s, z)\ d sdz 



p(s, z)h(s, z)d 2 sdz) 



(10) 



The well known approximation that is valid for intermediate and heavy nuclei for convolution 
of the profile function T s (b — s) with nuclear density [16J is 



T s (b — s)p(s, z)d 2 sdz 



a 



p{b, z)dz, 



(11) 



4tt 



where a' = — — / s (0). On general grounds it is known that the spin-nonflip nucleon amplitude 

ik 

in photoproduction of pseudoscalar mesons would vanish at zero angles. Thus, we choose 
for the elementary photoproduction amplitude on the nucleon the parametrization 



f P (q) = (n-q)(f(q). 



(12) 



k x e 



Here ip(0) ^ 0; n 

k 

section for photoproduction off the nucleon is — - = -g 2 |<^(g)| 2 , and the relevant profile 



with e the polarization vector of the photon. The relevant cross 
da 1 
~dn ~ 2* 



function becomes 



r P (b - s) 
r(b-s) 



-in 



dT{b - s) 
db 



2ixik 



^ b -^<p(q)d 2 q. 



(13) 



Accounting for the fact that the profile function T p (b — s) changes much slower than the 
nuclear density p(s, z), it is straight forward to obtain from Eqs. l9land[T3l 



h(s, z) 



_ .o< ^ dT(s,z) 

V$-*-2"'—9g- 



Lf(q)e 



iq-s+iAz 



x exp ( - tt^O, z] 

POO 

T(s,z) = I p(s,z')dz'. 

J z 



(14) 



Substituting this expression into Eq. [TO], we get for the incoherent cross section of the process 
under consideration: 



da.;. 



dtt 



ddo 
~dtt 



(q)[N(0,a) 



a 



p{s, z) 



A 

dT(s, z) 



ds 



x exp ( — crT(s, z)jd 2 sdz, 
F{q, A) = J e^ g+iAz p(s, z)d 2 sdz exp ( - ^T(s, z] 



TIG 



—J Ji(qs)p(s,z 2 

f (J 1 \ 

x exp y - —T(s,zi)J, 



9p(s,zi) 



ds 



sdsdz\dz2e 



iAzi 



N(0,a) 



1 — exp (—a f p(s, z)dz) 



d z s. 



(15) 



Let us mention the main^ifference between this expression and the one widely used for the 
incoherent cross section 



dcji. 



dn 



^(q)N(0,a) [1-G(t)\. 



(16) 



At small momentum transfer this expression goes to zero (G(0) = 1). On the other hand, this 
is not the case for Eq. [151 which differs from zero at q = 0. The complete suppression takes 
place only if one neglects the final state interaction of the produced meson (cr = 0). This 
fact has been well known for many years for incoherent cross section of elastic scattering 1 161 

nn 

and for diffractive photoproduction [171. 118I|. 



III. PAULI CORRELATIONS 



Up to now we have assumed the factorization for the nuclear ground state wave function 
(see Eq. [7]). As has been shown above, even in such simplified model the orthogonality 
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and completeness of final states nuclear wave functions (closure approximation) lead to 
suppression of the pseudoscalar mesons' yield at forward angles different from the common 
one due to the mesons' final state interactions in nuclei. 

Let us estimate the impact of the nucleon-nucleon correlations on the incoherent cross 
section. Here we are constrained by the correlations due to the Pauli exclusion principle, as 
their contribution is the largest one and they determine the behavior of the cross section at 
small momentum transfers. 

The ground state wave function of the nucleus can be cast in the form [3] 

p {A) (n,r 2 ...f A ) = |$i(fi,f 2 ...f A )| 2 = p(f l )p(f 2 )...p{f A ) 

+ Yl Wi)p(r 2 )...p(f A )] + (17) 

contraction 

where the pair contraction is defined by A(fi,r 2 ) = p( 2 '(r 1,7*2) — p(r\)p(r 2 ), and the single- 
particle and two-particle densities 

p{n)=J p [A \r 1 ,r 2 ...r A )d 3 r 2 d 3 r z ...d i r A , 

P {2) (r 1 ,r 2 ) = J p {A) (f 1 ,f 2 ...f A )d 3 r 3 d 3 r A ...d 3 r A . (18) 
Let us introduce the so called mixed density [3] 

A 

d(r 1 ,r 2 ) = J2<f>j(nW j (?2), (19) 

where <f)j{r) is the wave function of the single nucleon, so the nuclear density is given by 
p(r) = d(f,f). The two-body density can be expressed through the mixed density 

1 



p (2) (fi,r 2 ) 



A(A - T) K f i> f iM^2,r 2 ) 
d(Fi,r 2 )d(r ? 2 ,T ! i)]. (20) 



Substituting this expression into the general formula, Eq. [6], for the summed cross section 
and using the relation 

J d(fi,f 2 )d(r 2 ,f l )d 3 r 2 = d{n,r x ), (21) 



we obtain 



J" c = / d(s, z } s, z)\h(s, z)\ 2 d 2 sdz 



- J d(s 1 ,z 1 ,s 2 ,z 2 )d*(s 2 ,z 2 ,s 1 ,z 1 ) 

x h(si, Zi)h*(s 2 , z 2 )d 2 Sidzid 2 s 2 dz 2 . (22) 



It is straight forward 
the well known form 



;o check that in the Born approximation (er = 0) this expression takes 



G(q, A) = J e ^-ft)+*A(«i-«) 

X d(si,^i,S2,Z 2 )d*(s2,^2,Si,2;i) 

x d 2 Sid 2 S2dz\dz2. (23) 

To proceed with the calculations for the incoherent cross section, one has to choose the 
model for the mixed density. We will consider the photoproduction from light nuclei, for 
which the shell model with harmonic-oscillator wave functions works well. 

Let us first consider a closed shell nucleus, for instance 16 O. The mixed density in the 
independent particle model for this nucleus with harmonic-oscillator wave functions |7|] is 

fif*2 \ / r\ + r\ 



d ^>=^( i+2 ifjn-^r> (24) 

Substituting this equation into Eq. [22] we calculate the incoherent cross section. 

In Fig. [1] we plot the dependence of the incoherent cross section for photoproduction of 
7T° mesons off 16 nucleus as a function of the production angle. The oscillator parameter 
for 16 is ao = 1.81 fm {^(J, and the total cross section of ir°N interaction is o = 27 mb. 
The solid line is calculated by Eq. [[HI The long-dashed and dotted curves are the results of 
calculations according to Eqs. [15] and [22] respectively. 



For the 12 C nucleus, we choose the mixed density in the matrix form 211 ]: 



2 2 

d{fx, f 2 ) = „ . [l + 77-2 (2fi • f 2 + ia ■ (fi x r 2 )) 1 
agTr- 5 /^ L dag j 

xex P (-±^t). (25) 

In this case the single-particle nucleon density is expressed through mixed density by the 
relation 

p{r) = Tr(d(r,r)). (26) 
The results of our calculations for 12 C using mixed density, Eq. [251 with oscillator pa- 



rameter ao = 1.6 fm 16] are shown in Fig. [2] The solid line is calculated using Eq. [16] The 



long-dashed and the dotted curves have the same meaning as in Fig. [T] 
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FIG. 1: Incoherent differential cross section for 7r° photoproduction off 16 O nucleus as a function 
of production angle in the lab system. The long-dashed line is calculated by Eq. [15] with single- 
nucleon density relevant to harmonic-oscillator wave functions; the dotted line is the incoherent 
cross section taking into account Pauli correlations, Eqs. [22] and [Ml the solid line is the standard 
parametrization as in Eq. [TBI 

As seen from the figures, the correct accounting of pion absorption changes the behavior 
of the incoherent cross sections at forward angles which can be crucial for the extraction of 
radiative decay widths of pseudoscalar mesons. 

IV. SUMMARY AND CONCLUSIONS 

In recent years it has become possible to perform high precision measurements of photo- 
production cross sections for light pseudoscalar mesons on nuclei. The extraction of meson 
radiative decay widths with high accuracy from Primakoff type of experiments requires 
state-of-the-art theoretical descriptions of all participating processes. We have calculated 
the incoherent photoproduction cross sections for 7r° using the Glauber multiple scattering 
theory. On that way, we have revised the existing approaches for incoherent non- diffract ive 
photoproduction at forward angles and obtained new expressions (see Eq. IT5l) which cor- 
rectly account for meson absorption in the final state. This expression correctly takes into 
account the fact that the differential cross section on a nucleon is zero at zero production 
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FIG. 2: Incoherent differential cross section for 7r° photoproduction off 12 C nucleus as a function 
of the production angle in the lab system. The curves have the same meaning as in Fig. [TJ 

angle. Moreover, the zero in the nucleon cross section leads to additional terms in the in- 
coherent cross section without which the incoherent cross section would be negative at very 
small angles. 

Pauli correlations between nucleons on which photoproduction takes place are properly 
taken into account as their effects are dominant at forward angles. To account for Pauli 
suppression, the ground state wave functions of nuclei have to be antisymmetrized as the 
nucleons are fermions. This can be done if one works using the mixed densities, as in 

Taking into account the above considerations, the general expressions of the incoherent 
photoproduction cross sections for pseudoscalar mesons off nuclei were obtained for the first 
time (see Eq. I22l) . Using these new expressions and taking the mixed and single nucleon den- 
sities in correspondence with the harmonic-oscillator model, we have performed calculations 
for a closed shell light nucleus (oxygen) and an unclosed shell nucleus (carbon). 

V. ACKNOWLEDGMENTS 

This work was inspired by the PrimEx experiment at Jefferson Lab to measure the 7r° 
radiative decay width via the Primakoff effect with high precision. It was partially supported 
by the National Science Foundation grants PHY-0245407 and PHY-0555524. One of the 



11 



authors (S.G.) is grateful to Alexander Tarasov for useful comments. 



[1] H. Primakoff, Phys. Rev. 81, 899 (1951). 

[2] A. Gasparian, invited talk at the American Physical Society's (APS) April meeting, Jack- 
sonville, FL, April 2007; L. Gan, invited talk at the APS April meeting, St. Louis, MO, April 
2008. 

[3] S. Gevorkyan, A .Gasparian, L. Gan, I. Larin, M. Khandaker, arXiv:hep-ph/0903.4715 
(2009). 

[4] J. Goity, A. Bernstein, B. Holstein, Phys. Rev. D66, 076014 (2002). 
[5] B. Ananthanarayan and B. Moussallam, JHEP 05, 052 (2002). 
[6] B. L. Ioffe and A. G. Oganesian, Phys. Lett. B647, 389 (2007). 
[7] C. A. Engelbrecht, Phys. Rev. 133, 988 (1964). 
[8] G. Bellettini et al, Nuovo Cimento 66A, 243 (1970). 

[9] G. F. Chew, M. L. Goldberger, F. E. Low and Y. Nambu, Phys. Rev. 106, 1345 (1957). 
[10] A. Sibirtsev et al, arXiv : hep-ph/0902 . 1819 (2009). 
[11] G. Faldt, Nucl. Phys. B43, 591 (1972). 
[12] G. Faldt, Nucl. Phys. B62, 86 (1973). 
[13] G. Faldt and P. Osland, Nucl. Phys. B126, 221 (1977). 

[14] R. J. Glauber, in "High-energy Physics and Nuclear Structure", Proceedings of the II Inter- 
national Conference, Rehovoth, p. 311, 1967. 
[15] S. Gevorkyan, A. Tarasov, O. Zaimidoroga, preprint JINR P2-6581 (1972). 
[16] R. J. Glauber and G. Matthiae, Nucl. Phys. B21, 135 (1970). 
[17] E. Moniz and G. Nixon, Phys. Lett. B30, 393 (1969). 
[18] E. Moniz and G. Nixon, Annals of Physics 67, 58 (1971). 
[19] M. Braunschweig et al, Nucl. Phys. B51, 157 (1973). 

[20] R. Hofstadter and H. R. Collard in "Landolt-Bornstein" , New Series Group 1, vol. 2, Springer- 

Verlag, Berlin (1967). 
[21] A. Tarasov, Phys. of Elem. Part, and Atom. Nucl. (PEPAN) 7, 771 (1976). 



12 



